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Abstract 

We suggested a geometric approach to address the external field influence on the DNA molecules, described 
by the WLC model via geometric coupling. It consists in the introduction of the effective metrics depending 
on the potential of the external field, with further re-definition of the arc-length parameter and of the extrinsic 
curvatures of the DNA molecules. It yields the nontrivial impact of the external field in the internal energy of 
macromolecules. We give the Hamiltonian formulation of this model and perform its preliminary analysis in 
the redefinition of the initial energy density. 

1 Introduction 

In recent years, single-molecule techniques have evolved into a powerful toolset for studying the dynamical behavior 
of many biological processes. It is now possible to follow on the nanometer and millisecond scales the individual 
trajectory of a single enzyme as it catalyzes a reaction, a molecular motor as it translocates, or a single polypeptide 
or nucleic acid molecule as it unfolds and refolds [JJ [2] [3] . 

One of the widely used applications of single-molecule manipulations is the direct investigation of double- 
stranded DNA (dsDNA). The single- molecule stretching experiments gave a big contribution to the understanding 
of the dsDNA structure and functionality [4j [5] . The interpretation of any stretching experiment is closely related 
with the mechanism of elasticity. For example, the dsDNA stretching experiments are usually interpreted in terms 
of worm- like chains (WLC) [H [5J [3 H] , while the ssDNA or ssRNA are usually described by freely jointed chains 
with elastic bonds (EFJC) 0. 

The WLC Hamiltonian is presented as a functional on a differentiable space curve 7 of fixed length L, 



cWLC _ IpkpT 



/ Kfds , (1.1) 

J <v 



where r(s) is the radius vector of an arbitrary point of the curve as a function of the contour distance s from one 
end to that point, T is the temperature, fee is the Boltzmann constant, l p is the persistence length and «i is the 
first curvature (rigidity) of the curve, which is defined by the use of the unit vector ei(s), ef = 1 tangential to the 
curve: K\ = |dei/ds|, ei = where s is the arc-length of the curve r(r), 

dv 

\dr\, ds=\^\dr (1.2) 
dr 

An example of the spatial constraints, imposed on the polymer chain, is the helical structure observed in many 
biological and synthetic polymers [9]. In order to describe the helical structure of the chain space, following Bugle 
and Fujita [TU] (see also [IT]), we need to introduce the unit curvature vector e 2 = r^ygfj and define the potential 
energy of the chain as follows: 

£r = ^f j 'd*[(«i-a) a + (K 2 -£ 



(1.3) 



where k\ = e\ ■ e 2 is a first curvature (rigidity) of 7 and k 2 = y/ \de2 / ds\ 2 — k\ is a second curvature (torsion) 
of the curve 7. The a and j3 parameters are supposed to be the equilibrium bending and torsion constants. The 
absolute minimum of the potential energy is given by Ki = a and k 2 = j3 and corresponds to the regular helix 
with radius r = and pitch h = Jrj^l ■ 

In this respect the paper by Feoli, Nesterenko and Scarpetta [12] deserves to be mentioned. There the authors 
suggested to describe the DNA molecules by the functional 



£0= (co + ciKi)ds, (1.4) 
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whose minimum corresponds to the helical configuration of the curve. Unfortunately, these authors interpreted 
the functional (|2.1I) as a free energy of the DNA molecules, which yields some inconsistencies in the interpretation 
of the results. The importance of that paper is not only in the suggestion to describe the DNA molecules by the 
functionals linear on the curvature, but to operate with the functionals depending on extrinsic curvatures by the 
use of the Dirac's theory of constrained systems [13]. Further below, we will discuss this issue in more detail. 

It is not occasional, that in all listed models the energy of the macromolecule is defined by the use of extrinsic 
curvatures of the macromolecule curve. Indeed, the (linear) density of the energy of the macromolecule should a 
priori depend on the geometric characteristic of the curves and be independent on the way of their description. 
Hence, it should depend on the functions of the curve 7 which are invariant under its motions (rotations and 
translations), i.e. on its extrinsic curvatures ft/, < I < D — 1 (where D is the dimensionality of the space) [T3] 

From this viewpoint the common approach, for taking into account the external field influence, seems to be 
ill-defined. It consists in the naive adding to the energy functional given further below in (|2.ip of the potential 
energy term [15] 

£ = £0+ [ V(r)ds, (1.5) 

where V(r(s)) is the potential of the external field. 

The geometric ( "macroscopic" ) inconsistency is that the energy functional of the macromolecule becomes the 
function of a quantity, which is not invariant under the motions of the curve (rotations and translation in Euclidean 
space). 

The physical ("microscopic") inconsistency of this approach is the independence of the chain flexibility (e.g. the 
persistence length l p ) on the external field. While the stretching or the compressing fields substantially restricted 
the set of available conformations of each repeated unit of the polymeric chain. Thus, a strong enough external 
field should affect the chain flexibility, and the persistence length l p . This issue should be relevant, e.g. for the 
confinement-induced changes of the persistence length observed in |16j . The effective compressing field, describing 
the intra - chain attraction and chain confinement in restricted geometries yields a substantial change of the 
apparent persistence length pT] [18] . 

The purpose of the present paper is to suggest an alternative way for taking into account the external field 
influence, which seems to be free from the above discrepancies. Namely, we suggest to involve the external field 
in the definition of the line element ds, associating it with the effective metric 

g AB = (1 + V{v)f S AB = n 2 (r)S AB , (1.6) 

and to define the extrinsic curvatures hi, in accordance with Frenet equations, on the space equipped with this 
effective metric. 

Hence, we suggest to preserve the functional dependence of the energy density from the extrinsic curvature, 
but modify the definition of the arc-length, and consequently, of the extrinsic curvatures. 
The paper is organized as follows. 

In the Second section we present the explicit expression of the suggested geometric model of the macromolecule 
interaction with the external field, and discuss the analogies of this model with relativistic particle systems and 
with geometric optics. 

In the Third section we give its Hamiltonian formulation and analyze it in the framework of the Dirac's theory 
of constrained systems. 

2 The model 

We suggest to describe the external field influence on the macromolecule by the energy functional 

£ = / F(ki,... ., K N )ds, (2.1) 

J -y 

where 

• F(kj) is the energy density in the initial (noninteracting) DNA model. It can be chosen say, as in (jl.l[) , 

Q, (ED 
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• The element of arclength (or control length) ds is defined in accordance with the metric (ll.6[) . 



ds = n(r)\r\dr = sdr, (2.2) 
where r is an arbitrary parameter defining the curve 7: r = r(r) 
• The extrinsic curvatures kj are defined, via Frenet equations, on the space equipped by the effective metric 

da, 

dr De a 

T3-=ei, — — = n a e a+1 - K a -ie a -i. (2.3) 
sdr sdr 

Here 

-r- = ^-+f(t), (f)i=T^ c i c , e a -e b = 5ab, e = e D +i = 0, I = 1, . . . , D 
or dr 

with eje,' = e^ A gABe^ B , being £> the dimensionality of the space and Tg C the Christoffel symbols of the 
effective metric 

Tab = u aScb + n B 5cA - ticSab, ua = Oa logn. (2.4) 

Similarly to the Euclidean case, the curvatures «i, . . . , /t^-i are positive functions, and the highest curvature 
(torsion), Kd—i can take both positive and negative values. If some ki 7^ 0, then m ^ i = 1, 2, . . . , I — 1. 
Vice versa, from = 0, it follows that m = = I + 1, . . . , D — 1. 

Explicit expressions for external curvatures look as follows: 

v/dct.gz+idet (f r d J r . . 

K/ = — ^ — ' (5/) - = WW ^ =1 '---' J - (2 - 5) 

In particular, the first curvature reads 



Kl "^wY K (/^) i + ^f^' (2 - 6) 

where Kifi a t)x is the first curvature of the curve on the flat space. Hence, within the suggested approach, we have 
got the geometrically consistent energy functional of the macromolecule, where the influence of the external field 
in the internal energy is taken into account. Since we deal with the macroscopic model (or a phenomenological 
description) of the macromolecule, we can assume that n(r) parametrically depends on the temperature T. For 
the correspondence with the existing models we can choose, e.g. n — n(T,r) = fcsT(l + ^f)- The environment 
influence on the chain elasticity and flexibility is usually described in terms of the persistence length. For example, 
to describe the DNA behavior in nanochannels and nanoslits, the concept of the global persistence length has 
been introduced in [19] , as a typical distance between hairpins in the confined chain. The main advantage of 
the obtained result is the deep non-trivial relationship between the bending elasticity and the external field. The 
proposed approach seems to be useful for a better understanding of the mechanism of elasticity of macromolecules, 
in stretching and confining fields. 



The price for the geometrization of the model is the further mathematical complication of the system. In spite 
of the apparent simplicity, the analysis of models with the energy functional (|2.1|) is not a simple task: the matter 
is that the energy density F in (|2.1|) depends on higher derivatives and should be analyzed in the framework of 
the Dirac's theory of constrained systems |13j . Fortunately, in the eighties such functionals were widely studied 
in the context of the (relativistic) spinning particles. Namely, the functionals (|2.1[) in the Minkowski space- 
time, with s being a world-line parameter were considered. These studies were inspired by the papers Polyakov 
devoted to the Chern-Simons theories and rigid strings [20]. Various systems depending on the first and the 
second curvatures of the path in three- and four- dimensional Minkowski spaces were studied in details, mostly 
by Plyushchay (see (21] [22] and references therein). By this reason the relevant formalism for analyzing such 
systems on the base of the Frenet equations has been developed by one of the authors [23] . Almost all systems 
depending on extrinsic curvatures either by isospin (i.e. by internal degrees of freedom), or spin (in the case when 
they are proportional to a single curvature, see [22l [24]. In the context of macromolecules it means, that the 
systems with energy density depending on extrinsic curvatures of macromolecules should have internal degrees 
of freedom. Now, any macromolecule is a sequence of atoms connected by chemical bonds. That is why, the 
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universal way of any change of chain conformation is the change of dihedral angles, caused by rotation around 
the chemical bonds. Thus, any polymeric chain has its internal rotational degrees of freedom, which should be 
obviously identified with the isospin of the system. These studies were performed on flat spaces, mostly. The only 
exceptions are, seemingly |25[ 124]. A qualitative observation there is that, for Lagrangians (in our case: energy 
densities) with linear dependence on extrinsic curvatures, the introduction of non-Euclidean metrics drastically 
changes the properties of the system, including the dimensionality of the internal space. Only for the constant 
curvature spaces (e.g. spheres, hyperboloids) non-Euclidean metrics do not change essentially the properties of 
such systems. In three-dimensional space, in our notation it corresponds to the choice 

n(r) =a /(l + /3r 2 ) (2.7) 

corresponding to the three-dimensional sphere or hyperboloid (depending on the choice of the sign of /?) . 

The suggested DNA model has an obvious analogy with the geometric optic, where the trajectory of light in 
the media with refraction index n(r) corresponds to the minimum of the functional (under assumption that the 
helicity of light is neglected) 

SFermat = Jds, ds = n(r) |r|d<r. (2.8) 

Hence, the generic action (|2.ip defined on such "effectively" curved space, will take into account the interaction of 
the spin (helicity) of the light with the media, i.e. the phenomenon of the "optical Hall effect" , consisting in the 
deviation of the light trajectory from that given by Fermat Principle, due to the feedback from the polarization 
[2"5] . Let us also mention, in this respect, the recent paper [37], where the analysis of the model (|1.4I) with a 
non-constant effective metric has been performed in some details, in the context of geometric optic. 

So, we suggest to describe the DNA molecules interacting with the external field by the functional (|2.ip . where 
the arc-length is defined by the expression ds = (1 + V(r))|r \dr, and the extrinsic curvatures Kx, . . . , kn, N < D 
are defined in accordance with modified Frenet equations (|2.3|) . 



3 Hamiltonian formulation 

It is seen from the Frenet equations (|2 . 3[) . that the I-th extrinsic curvature ki depends on the (I + l)-th derivative, 
= ftr(rw)), so that the energy density F(Ki) in (|2.ip depends on the (N + l)-th derivative. By this reason 
the study of its Euler-Lagrange equations is not a simple task. To simplify the analysis of (|2.ip . let us consider 
the variationally equivalent functional, whose integrand depends on the first derivatives only (see [23j GH [28] for 
further details). 

Precisely, we replace the initial energy density C by the following variationally equivalent one: 



F = F(ki, . . . , kn) - A \ kn 



>(De 



N 



V sdr 



"N-l 



+ P ( z— - e i ) + y^P'- 1 1 1 _ K «-i e ' + K t - 2 <ii-2 J - V"d lJ (e^- - 6ij) , (3.1) 
\sdr J r-J V sdr J —. 

where s, d lJ , Pi-i, ej, A play the role of independent variables, and the relation 

2 _ ( De N \ 2 2 

K N — I gfjlj- ) ~ K iV-l' 

which follows from (|2.3p is taken into account. Here and further below, we use the following notation for the scalar 
products: pt ■ Pj = P(i)A9 AB P(j)B: Pi ' e j = P(i)A e fi where, upon summation the indices A, B take the values 
1,...,D, and i,j = 1,...,N. 

It is easy to see, that varying the energy density F by kn, we shall get A = dF/dn^, while its variations by 
p, Pi_i and dij will yield the Frenet equations restoring the dependence of Ki on r*- 1 '. 

The Hamiltonian formulation of the system with the singular "Lagrangian" (|3.ip on the Euclidean space was 
carried out in [23 . Its extension to curved space is straightforward. For the "Lagrangians" linear on extrinsic 
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curvatures it was considered in 23], and for the generic case in [35]. In accordance with this paper, the Eulcr- 
Lagrange equations of the functional with the energy density F are given by the constrained Hamiltonian system, 

jv 

lun = dp A dr + dpi A de.;, H = ho + 2_. K i<t>i + ^ & ( e i e j ~~ (3-2) 

i—l i i.j 

with the following set of constraints (including some gauge fixing conditions, see for details [2"31 12"5] ): 

N 

00 = irei + K i F ,i ~ F ~ °> where * = p - I\ r A = ^ T ABP{i)C&f (3.3) 

i i—l 

Piej w 0, i > j (3.4) 
- w 0, (3.5) 

= Pi-iej - p,ej_i - F,:_i w 0, 2(j> N = -=— [ PtvPat - y^(pAre i )(p A re i ) ] - Fjv « (3.6) 



For det d 2 F/dKidKj 0, we can resolve the equations <\3.6\i , expressing m via Pi,ej. After that we should 
substitute these expressions for Kj in p.3[) . 

Otherwise, for corank d 2 F '/ 'dKidKj = M ^ 0, we can resolve only (TV — M) equations in (|3.6p . while the 
rest will appear as a primary constraint, with the M undefined functions on K{ playing the role of Lagrangian 
multipliers. 

After these manipulations we will get the constrained Hamiltonian system with N(N— 1)+1 primary constraints 
given by (|3.3|) - (|3.5|) and with M primary constraints obtained from (|3.6p . Then, we should perform the Dirac's 
procedure of "stabilization" of the above set of constraints, i.e. require {. . . , H} « 0, where by ... we encoded 
the set of primary constraints. This requirement will either fix the values of Lagrangian multipliers, or yield an 
additional set of "secondary" constraints. Then, stabilizing the obtained set of secondary constraints, we shall 
either get the new set of secondary constraints, or will further fix the values of Lagrangian multipliers, and so on 
(see, e.g., [13]). 

It is easy to observe, that the "stabilization" of the constraints (|3.4j) fixes the values of the Lagrangian multipliers 
d lJ , while the stabilization of the remaining set depends on the concrete properties of the function F(kx, . . . , kn)- 

For example, for the energy densities with a non-degenerate quadratic dependence on curvatures, F — cijKiKj , 
detcy ^ 0, which involves, as a particular case, the WLC model and Buggle-Fujita model (| 1 . 3[) . there are no 
secondary constraints. The whole set of (primary) constraints is divided in the second-class constraints given by 
(|3.4p . (|3.5[) . and the single first-class constraint given by (|3.3j) . As a consequence, we get that the dimensionality 
of the physical phase space is 2D(N + 1) - 2 - N(N - 1) = 2DN - N(N - 1). In contrast with this case, in 
the systems with linear energy densities, F — CiKi + Co, the whole set of constraints and their algebra, essentially 
depend on the values of the constants Cj and on the properties of the effective metric (11.61) . 

Hence, in contrast with the functionals with linear dependence on curvature, in the case of quadratic depen- 
dence, the introduction of the curved metric does not essentially change the structure of the phase space of the 
system, being reflected only on the functional dependence of the Hamiltonian from the external field. 

To complete this general consideration, let us present the following useful expressions, which should be useful 
for the further analysis of the systems under consideration. 

At first, let us write down the non-zero Poisson brackets between the functions 7re,;, 7rpj, 7r 2 : 



{7rei,7rp } =TZ(e i ,p :j ) --K 2 5ij, {-Tre^Tr 2 } = -27^(7r,ei), {ttp^tt 2 } = -21l(-K,Pi) 



2\ _ _o<r>/v ™^ (3- 7 ) 



where 

N 

7l(a,b) = 5^ii(p i |e 4 ,a,b). (3.8) 

j=i 

is the Riemann tensor associated with the effective metric (|1.6p . 

Then, let us present the expression for the Riemann tensor in the three-dimensional space, where it is uniquely 
expressed via the Ricci tensor 

Rq 

Rabcd = R(2)ac9bd — R(2)ad9bc + R(2)bd9ac — R(2)bc9ad + (9ad9bc — 9ac9bd) , (3.9) 
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where i?2 is the Ricci tensor and i?o is the scalar curvature of the space. For the conformal fiat metric (ll.6[) the 
Ricci tensor is defined by the formulae 



Rac = ~9Ac{nDD + 3n|,) - n AiC + n A n C) n A = d A \ogn. (3.10) 

All considerations above were quite general, with a generic functional dependence of the energy density on the 
extrinsic curvatures, and they have been formulated in the D-dimensional space. As an example, let us give the 
explicit formulation of the WLC model (jl.ip . with the additional constant term Cq 

F=^kI+co. (3.11) 

In this case the Hamiltonian system (|3.2p reads 

uj = dp Adr + dpi A de%, H = nei + — — cq + ( c — -7rei ) (ei • ei — 1) , (3-12) 

2c \ 2 J 

with the following set of constraints: 

P 2 

ei • ei - 1 « 0, pi • ei w 0, 7rei H — c « 0. (3.13) 

2c 

The extrinsic curvatures looks as follows: 

[rrf If 3 \ 

«i = V^' K 2 = ^"3 7T 2 -4c 7rei + -(Trei) 2 ) ~ k\. (3.14) 

V C C \ 4 J 

The evolution equations for these curvatures could be straightforwardly obtained by the use of Poisson bracket 
relations (|3.7p . They will obviously depend from the Riemann tensor p. 81) . 

We suppose to present the detailed study of the WLC model in a forthcoming paper. 



4 Discussion 

The statistical mechanics of the WLC has been developed by regarding it as a differentiable space curve of fixed 
length [H],[TT]. If r(s) is the radius vector of an arbitrary point of the curve as a function of the contour distance 
s from one end to that point, persistence length l p can be defined as a correlation length of the tangent vectors 
fluctuations [9] 

(u(«)u(« + r)) ~ e- r ^ lp , (4.15) 

where u(s) is the unit vector tangential to the curve at the point s. Thus, l p is statistical quantity strictly dependent 
on the set of available chain conformations. It is well known that the intra-chain interactions, e.g. electrostatic 
repulsion |29j or geometrical restrictions imposed on the chain conformation |17j renormalise persistent length l p . 
However, these particular issues addressed before could not provide the systematic approach to the influence of 
interactions and external fields on the flexibility and elasticity of macromolecules. 

In our opinion just the artificial separation of the terms describing chain elasticity and external fields or 
interactions presented in (|1.5|) results to the strong necessity of the persistence length renormalization. Besides 
of artificial renormalization of the persistence length the commonly used additive scheme produced numerous 
inconsistencies with the description of chain flexibility. For example, stretching experiments with dsDNA of length 
about 0.6 — 7 fim produced the resulting fitted value for the persistence length depended on the contour length of 
the dsDNA [30] . The elastic energy of highly bent dsDNA conformations is lower than predicted by WLC model 
[31] . The list of inconsistencies can be continued. 

The proposed above approach permits to take into account the external fields (or interactions) effect on the 
chain rigidity and elasticity directly without to invoke the artificial renormalization of persistence length l p and 
promised to extend our understanding of the physical properties of biological macromolecules. It consists in the 
introduction of the effective metrics depending on the potential of the external field, with further re-definition 
of the arc-length parameter and of the extrinsic curvatures of the chain. These redefinitions obviously can be 
interpreted as an impact of the external field in the internal energy of macromolecules, and in the natural changes 
of the persistent length. The proposed way to address the effect of the external fields on the conformational 
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statistics of macromolecules seems to be much more general than commonly used. Surely, we do not get any 
concrete prediction on the physical relevance of the suggested coupling. But we believe in its correctness, since it 
was suggested in a first principles in accordance with quite general symmetry considerations, and from this point 
of view, can not be wrong. While common scheme of counting external field influence coupling is pure heuristic 
one, and has no any motivation, except formal simplicity (or analogy with nonrelativistic mechanical systems). 
Moreover, any other coupling scheme, based on the "microscopic consideration" , will be in our opinion a priory 
more rough and simplified, than presented one. The proposed symmetry - based approach to the external field 
influence is in the same spirit as Bugl and Fujita model |10j taking into account bending and torsion in terms of 
the first and second curvatures. 

Technical attractivity of the suggested DNA model insist in the with spinning particles on curved space, as 
well as with geometric optics. This allows us easily to evaluate qualitative properties of the systems, as well as to 
involve, in the study of the DNA macromolecules, the tools of symplectic and Riemann geometry, Dirac's theory 
of constrained systems and of the rclativistic quantum mechanics. 
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